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We examine antiferromagnetic and d-wave superfluid phases of cold fermionic atoms witii repulsive 
interactions in a two-dimensional optical lattice combined with a harmonic trapping potential. For 
experimentally realistic parameters, the trapping potential leads to the coexistence of magnetic and 
superfluid ordered phases with the normal phase. We study the intriguing shell structures arising 
from the competition between the magnetic and superfluid order as a function of the filling fraction. 
In certain cases, antiferromagnetism induces superfluidity by charge redistributions. We furthermore 
demonstrate how these shell structures can be detected as distinct antibunching dips and pairing 
peaks in the density-density correlation function probed in expansion experiments. 



PACS numbers: 03.75.Ss, 05.30.Fk, 74.25.Ha, 74.72.-h 

The trapping of ultracold atoms in optical lattices 
opens up the possibility to study quantum systems in 
periodic potentials with unprecedented experimental ver- 
satility. One can mimic strongly correlated systems rele- 
vant for paradigmatic condensed matter applications as 
well as creating entirely new structures. The pace of ex- 
perimental progress is impressive. The quantum phase 
transition between a superfluid and a Mott insulator has 
been observed for bosons [H, and Fermi surface effects 
were reported for fermions 2] . Bunching and antibunch- 
ing effects in the density-density correlations were found 
for bosons and fermions in expansion experiments [3|1 and 
evidence of s-wave pairing was recently presented 

A major goal is to study two component [a = ±) 
fermionic atoms in a two-dimensional (2D) optical lat- 
tice with repulsive interactions. Such a system is well 
described by the Hubbard model whose phase diagram is 
controversial and directly related to the physics of high- 
Tc superconductors. For filling fractions close to one par- 
ticle per site (n = 1), the system is antiferromagnetic 
(AF) whereas for smaller densities the true ground state 
is unknown. An important question is whether a d-wave 
superconducting state arises solely from repulsive inter- 
actions, but many other open questions remain, including 
suggested stripe and checkerboard charge ordered ground 
states jsj that possibly coexist with superconductivity for 
some range of doping x = 1 — n. Experimentally, there 
is strong evidence for such inhomogeneous states in the 
cuprates0, @, 0|, but it is unclear if they are important 
for pairing. An aim of the cold gas experiments is to 
improve our understanding of this complicated problem. 

Motivated by this, we study a two component gas of 
fermionic atoms in a 2D optical lattice. We include an 
external harmonic potential which originates from the 
Gaussian profile of the laser beams generating the trap. 
A main purpose of this paper is to study the interplay 
between this harmonic trapping potential and two of 
the most dominating ordered phases of the homogeneous 
Hubbard model: the antiferromagnetic and d-wave su- 



perconducting phases. Previous theoretical studies have 
focused on the one-dimensional (ID) case [1]. For real- 
istic parameters, we find that the magnetic and d-wave 
superfluid (dSF) phases coexist and form shell structures, 
in analogy with what has been observed for bosons 0]. 
We then examine how these shell structures can be de- 
tected in density-density correlations by expansion ex- 
periments similar to those performed for ideal gases Q . 

A two-component Fermi gas in an o ptic al lattice is well 
described by the Hubbard model [lol. Illj 



1 



-mcu' 
2 



(1) 



where di^ are annihilation operators for localized atoms 
on site i [at position = {Xi,Yi)] with spin a and 
fiia = fllo-flio- is the number operator. The parameters 
t and U > denote hopping between nearest-neighbor 
sites (ij) and on-site repulsion, respectively. They can 
be obtained from the atom-atom scattering length and 
the lowest band Wannier state 



The last term in 
Eq.iU]) describes the harmonic trapping potential. We 
use unrestricted Hartree-Fock approximation to decouple 
the interaction term in Eq.([T]). 

In addition to the AF state at half-filling, numeri- 
cal methods indicate that the repulsive Hubbard model 
prefers d-wave superconducting order in a certain region 
of the U — X phase diagram[13j. Therefore, in a har- 
monic trap we expect that AF and superfluid order may 
coexist in the atomic cloud. To model this we include 
explicitly a BCS d-wave term so that the final mean-field 
Hamiltonian becomes 

H ^-t a\/aj„ - ^ ^ + ^mcj^ ^ R\h,„ (2) 
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It is well-known that a nearest-neighbor attraction can 
generate such a d-wave BCS term at the self-consistent 
level IJ] . We stress that even though our mean-field 
model Eq.([2) is phcnomenological, it captures the ex- 
istence and competition of ordered phases at T = such 
as AF and superfluidity, and has been widely used pre- 
viously in the high-Tc community Note that even 
in ID it is possible to obtain qualitative information 
from mean-field theory, such as e.g. the shape of the 
density profiles and the existence of antiferromagnetic 
correlations 0. This gives further confidence in the 2D 
T = mean-field results presented below. 

Equation ([2|) can be diagonalized by the transforma- 
tion, a]^ = 'EE„,>oi<Ai + crvna{-i)^n-cr] yielding 
the Bogoliubov-de Gennes equations 
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The diagonal blocks are given by /C; 
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U{ni^a))Sij, where Vi — ^niuj^R^, and Sij and (5^^) 
are the Kronecker delta symbols connecting on-site and 
nearest neighbor sites, respectively. The off-diagonal 
block is Vij = 5(^ij^Aij, where in the homogeneous case 
Aij — -f(— )A on the a;(j/)-links corresponding to bulk 
da;2_j^2-wave pairing symmetry. Below, we restrict the 
discussion to T = and enforce self-consistency through 
iteration of the relations (fiia) = X^n kno^lOP ^-^^ ^ij — 

V'd(a^Tail-a^i%T) = E„KT(^Xi(J)-w«T(jXi(«)]■ 
Here, Vd is a coupling constant which, in principle, is a 

function of U, but at the phenomenological level becomes 

an independent parameter. 

We now present our results of the numerical solution of 
the mean-field equations varying the number of particles 
N trapped in a potential with ^muj'^cP /t = 0.025 where 
d = 1 is the lattice spacing. This yields experimentally 
realistic lattice sizes of the order ^ 40 x 40. For these 
relatively small systems and the parameters used in this 
paper, the superfluid coherence length ^ is comparable 
to the characteristic length scale A of the ordered phases. 
Therefore, we include the trapping potential exactly since 
the local density approximation cannot be expected to be 
valid because it assumes that ^ <C A. 

We first discuss the situation with no dSF order, i.e. 
Vd = 0, and use an x square lattice (A = 44) with 
open boundary conditions. In Fig. [T]we plot the density 
profile rii = {ni^) + {nii) and the staggered magnetization 
TO, = - (nj^))/2 with U/t = 4.0 for a 

varying number of trapped particles. The dashed lines 
in the left column display the density profile for an ideal 
gas {U = 0). For sufficiently high density (top panel), the 
center region is a band insulator in both the interacting 
and non- interacting limits. We see that the main effect of 
the interaction is the formation of AF regions at densities 
n ~ 1. Since the system is inhomogeneous due to the 
trapping potential, the AF order coexists with the normal 
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FIG. 1: (Color online) Left: line cut of the density (blue, solid 
lines) and trap potential (red, dot-dashed lines) through the 
center of the (spherically symmetric) trap, U/t = 4.0. The 
black dashed lines show the density profile for the noninter- 
acting case, (7 = 0. Right: the associated real-space staggered 
magnetization. The total number of fermions in the 44 x 44 
system are (top to bottom): 1936, 968, 484, 242. 



phase. This leads to steps in the density profile as the 
magnetic correlations favor n ~ 1. At the same time, 
the density is reduced in the center of the trap and the 
atomic cloud becomes more extended. Upon reducing 
the number of particles in the trap, the magnetization is 
seen to evolve from a ring structure to a center island. 

For larger systems, spin-density waves with ordering 
vectors other than the conventional AF Q — (tt, tt), e.g. 
stripe phases known from high-Tc materials, may become 
evident at filling away from n ~ 1. In that case we ex- 
pect the same overall results as those in Fig. [T]but with 
magnetism existing for a wider doping range. 

Next we discuss the possibility of superfiuid order fo- 
cusing on both the spatial distribution and the interplay 
with the magnetic order. In Fig. [2 we show the d-wave 
order parameter (left) defined as = (l/4)(Ai_i+£ -I- 



A,: 



A 



A 



where x{y) are the unit vec- 



tors along the x{y) axis, and the magnetization (right 
column). The density profiles (not shown) are similar to 
those presented in Fig. [TJ The spatial structures depicted 
in Figl2] come from the interplay between the trapping 
potential and the magnetic and superfluid correlations 
and can be understood as follows. The amplitude of A 



3 




0.2 



0.2 



I 



20 40 
40|||liiiiilllliil'l 

JlliF'liilsl 

20 li[jSL JSg gj 







jiiini 
































'is 










is 








■■■■■■■■MIHMI 
































■■ 














1 


ii 








ni 








. M 


■ 


m 


PIPP 


IB 


1 


1 


== 




1 


IHIIIII 


1 








20 40 




FIG. 2: (Color online) Left: spatial dependence of the d-wave 
superfluid order Ai with Vd/t = 2.0 and U/t = 4.0. Right: 
staggered magnetization obtained for the same set of parame- 
ters. The total number of fermions are (top to bottom): 726, 
484, 242, 144. 



peaks at one particle per site just like the AF order. This 
leads to a competition between dSF and AF order in re- 
gions around n ~ 1. We have chosen V^/i = 2.0 giving 
A ~ 0.1 5t (and a coherence length of a few lattice spac- 



ings) consistent with the numerical results of Ref. [13. 
For the resulting ratio VdjV — l/'2,, antiferromagnetism 
dominates and dSF order is generally left to exist in re- 
gions surrounding the AF order. This is seen, for exam- 
ple, in the two middle rows of Fig. [2] where AF exists in 
a center island with n ~ 1 surrounded by a ring of dSF. 
This is also the origin of the outer rim of the remarkable 
double-ring structure shown in the top left subplot of Fig. 
[2j There, however, the center dSF island is a case where 
magnetism surprisingly promotes dSF because of charge 
redistributions: The AF order decreases the density in 
the center of the trap below the threshold for generating 
dSF. For sufficiently small fillings, only a superfluid cloud 
is left in the trap (bottom row in Fig. [2]). 

We now address the important question of how the 
shell structures presented in Figs. [T] and [2] can be de- 
tected experimentally. Since magnetic and superfluid or- 
der coexist with the normal phase, it is a priori unclear 
how strong their experimental signatures will be. A well 
known experimental technique is to measure the density- 
density correlations of an expanding gas after the lattice 



has been switched off. The density-density correlation 
function (n(r)n(r')) after expansion time t probes the 
momentum correlation function (riqnq' } for the interact- 
ing system before the expansion with riq = a^^aqcr 
and r — qt/m [3]. We focus on the correlation function 

C(q,q') = ('^q'^q') - (?lq)(»T-q') = ^q,q'(?lq) (4) 

~l ^q'ng-^q"a'l + ^q'na'^qna^'i'na-bq'ma , 

na nrna 

where aqncribqna) = (1/iV) Wi„cr (wmo- ) exp(-iq • ) . 
For non-interacting fermions, C(q, q') has antibunching 
dips for q — q' = {nx^n , ny2T:) where and Uy are in- 
tegers as was recently observed f^]. Atoms forming an 
AF state will in addition exhibit antibunching dips in 
C(q, q') for q — q' = (jIxTt, riyT:) with odd integers n^^ riy, 
reflecting the period doubling due to the magnetic or- 
der. For a system without an external trapping potential 
characterized by an AF order parameter m — \ {riii ) — 
{nii)\/2, mean-field theory yields {riqUqi) ~ (f^q}("-q'} — 
E^(«q<Taq'^)(aq/^aq^) with (aj^^ aqv) = {Um/ 2)/2Eq 
for q = q' + (tt, tt) and Eg = y^e^ + {Um/2y. Here, 

eq = —2t[cos{qx) + cos{qy)] is the usual tight-binding 
spectrum. Likewise, for a homogeneous superfluid d- 
wave state, BCS theory yields (nqrtq') = (rtq)(rtq') + 

Sct l(aqTaq'-cr)P with |( aqa.a-q-cr)| ^ \Aq\/2Eq, for 

q — — q' and Eq = ^ (e^ — fi)^ + where Ag — 

2A[cos(g^) — cos(qj^)] is the d-wave gap. Pair correla- 
tions with s-wave symmetry were measured on the Bose- 
Einstein condensate side of a homogeneous system with 
a Feshbach resonance [l6l |. 

In order to examine how the antibunching dips and the 
pairing correlation peaks show up for the trapped lattice, 
we calculate C(q, q') for the shell structures shown in Fig. 
[T]and[l In Fig. [S^top left), we show a 2D map of C(q, q') 
with q' — {tt/2, tt/2) for the same parameters used in Fig. 
[T] (third row). Figure [Sjtop right) shows a cut along the 
diagonal qx—q'^ = qy — q'y Here, in addition to the usual 
lattice dipsjl], we clearly see the signature of the AF state 
in the additional antibunching dips at qx—Qx = (i""! ±7r). 
The momentum q' = (7r/2,7r/2) is close to the Fermi 
surface and hence maximizes the ratio of the AF dips to 
the 2tt lattice dips. Note that the perfect periodicity, i.e. 
the equal amplitude at points q ^ q + {nx2Tr , ny2Tr) , is 
an artifact of the model which neglects the spatial extend 
of the lowest Wannier function. 

We next discuss C(q, q') with both AF and dSF order 
present, and focus on the same particle filling as above 
which corresponds to the second row in Fig. O We show 
in Fig. [Sl^bottom left) momentum cuts in C(q, q') along 
a diagonal line with qx—q'x = Qy — Qy and q' = (7r/2, 7r/2) 
(solid line) and along a horizontal line qx — q'x with qy = 
and q' = (7r/2,0) (dashed line). Here, in addition to 
the AF dips, it is evident that the dSF phase displays 
bunching at q = — q' since the peaks appear at qx—q'x = 
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FIG. 3: (Color online) (top row) Density-density correlation 
function C versus q — q' (left) with q' — (7r/2,7r/2), and a 
cut (right) along the diagonal Qx — q'x ~ Qy — q'y Parame- 
ters are identical to those used in Fig. [1] (third row), (lower 
left) C(q, q') along a diagonal with Qx — q'x = Qy — q'y and 
q' = (7r/2, n/2) (solid line) and along a horizontal cut Qx — q'x 
with qy — and q' = (7r/2,0) (dashed line) for the same 
parameters as in Fig. [2] (second row), (lower right) C(q, q') 
versus q when q' — — q. 



Qx — tt/2 — — TT and ±27r displacements thereof. In Fig. 
[3]Jbottoni right) we fix q = — q' and plot a 2D map of 
C(q, — q). This clearly reveals the c?-wave symmetry of 
the pairing: the bunching is maximum along the x and y 
directions and minimum along x = iy. Also, the pairing 
is maximum around the Fermi surface as expected. 

The dSF pairing peaks in Fig. [3] are small and may 
be easier to detect at lower filling fractions where AF 
order is absent and a larger fraction of the particles par- 
ticipate in the pairing. Note that the AF antibunch- 
ing dips are also present in Fig. [3lJbottom right) at 
q = — q' = (±7r/2, ±7r/2). This is an example of how 
the co-existence of the magnetic and superfluid order 
can be detected as the presence of both AF antibunching 
and dSF bunching peaks in expansion experiments. The 
AF dips disappear for lower fillings when the Fermi sur- 
face moves below (±7r/2, ±7r/2) and the density vanishes 
there. For the same reason there are no lattice dips at 
(±7r, ±7r) in this image. 

In summary, we have studied the magnetic and super- 
fluid phases of repulsive fermionic atoms on a 2D lat- 
tice combined with a harmonic trap for experimentally 
realistic parameters. The Hubbard correlations result 
in intriguing magnetic real-space shell structures which 
may coexist and compete with a superfluid phase. These 
phases show up as distinct antiferromagnetic antibunch- 
ing dips and superfluid bunching peaks with d-wave sym- 



metry in the density-density correlations which can be 
probed in expansion experiments. The results are rele- 
vant to current experiments on atoms in optical lattices. 

B. M. A. is supported by the Danish Technical Re- 
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